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ABSTRACT 

We present optical measurements of the faint end of the luminosity function in the core of the Coma 
cluster. Dwarf galaxies are detected down to a limiting magnitude of R ~ 25.75 in images taken with 
the Hubble Space Telescope. This represents the faintest determination of the Coma luminosity 
function to date. With the assumption that errors due to cosmic variance are small, evidence is found 
for a steep faint end slope with a < — 2. Such a value is expected in theories in which reionization 
and other feedback processes are dependent on density. 

Subject headings: cosmology: observations — galaxies: clusters: individual (Coma) — galaxies: dwarf 
— galaxies: formation — galaxies: individual (NGC 4874) — galaxies: luminosity 
function, mass function 



1. INTRODUCTION 

The galaxy luminosity function (LF) is a fundamen- 
tal observational probe of galaxy formation and evolu- 
tion. Defined as the number of galaxies per unit magni- 
tude per unit area on the sky 2 , the LF depends on both 
the initial density fluctuation spectrum and on baryonic 
processes such as cooling, star formation and supernova 
feedback. Any theory attempting to explain how galax- 
ies form and evolve must test its predictions against the 
observed shape of the LF. 

The faint end of the LF is of particular interest. 
The standard cold dark matter (CDM) model predicts 
a steep faint end slope to the galaxy LF; the Press- 
Schechter approximation (1974) combined with a CDM- 
like power spectrum leads to an incr easing number 
of da rk halos as mass scales decrease (|White fc Reesl 
1978). A steep faint end slope has been seen in 
some environments, mainly those of high density (e.g., 
lTrenthaml[l998at IPhillips et all 11998(1 . In less extreme 
environments, however, fla tter sl opes are observed (e.g ., 
iPritchet fc van den Berghl 119991 iTrentham et al.ll200lj) . 
Why are CDM-consistent slopes not seen in all environ- 

1 Current address: Centre of the Universe, National Research 
Council of Canada, 5071 West Saanich Road, Victoria, BC V9E 
2E7, Canada; Margaret. L.Milne@nrc-cnrc.gc.ca 

2 More precisely, this is the definition of the luminosity distri- 
bution. The luminosity function is the number of galaxies per unit 
magnitude per unit volume. The term "luminosity function" is 
often used for both quantities in the literature. 



ments? 

Studies of the faint end of the LF are hampered by the 
extreme difficulty of obtaining complete samples of faint, 
low surface brightness galaxies. It is always possible that 
the differences in the faint end slope of the LF are simply 
due to faint galaxies being missed in less de nse regions. 
Other explanations have also been suggested. iTullv et all 
(2002), for example, use semi-analytic models to show 
that more low mass halos form earlier in regions that 
eventually become massive clusters. They then point 
out that the reionization of the universe could inhibit gas 
collapse in low mass halos. Combining these ideas, they 
conclude that dense regions start forming small, faint 
galaxies before reionization, and so can form many, while 
less dense regions begin small halo formation later, and 
so can only form few. 

In this work, we examine the faint end of the LF in the 
core of the Coma Cluster. Coma is a very rich cluster 
(Abell class 2); Tully et al.'s model would therefore pre- 
dict Coma to have formed many low mass halos before 
reionization. We expect to find a steep faint end slope to 
the Coma LF. We choose to use the me thod of statistical 
background subtraction l|Zwickvl|1957|) to remove back- 
ground galaxy counts from our sample of cluster galaxy 
counts: we count the number of galaxies in the cluster 
direction and then subtract the count of galaxies in a 
"blank sky" control field. 

Many Coma LFs have appeared in the literature (Ta- 
ble ^| . This work differs from previous studies in that 
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we use HST images of Coma and of our control field to 
determine the LF. The increased depth and resolution of 
these images improves the determination of the faint end 
slope of the Coma LF. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Coma Data 

Our Coma data consist of observations obtained from 
the HST archive. We selected 16 F606W exposures, to- 
taling 20400 s, and 6 F814W exposures, totaling 7800 
s, of a single WFPC2 pointing taken with the center of 
the PC chip placed on the nucleus of NGC 4874. Due 
to light contamination from this galaxy, we discarded 
the PC chip from this project. The combined area of 
the three WF chips is ~ 16650 arcsec 2 (after trimming), 
with a plate scale of 0.1" pix" 1 . For further details of 
these observations, see iKavelaars et aTl ll2000j) . 

The raw data were processed with the standard HST 
pre-processing pipeline. We then registered the images 
to within ~ 0.2 pix using ~ 20 stellar objects on each 
chip. We scaled the images based on their exposure time 
and applied a standard cosmic ray rejection algorithm. 
Finally, we combined the images using the median pixel 
value, ending with coadded images with effective expo- 
sure times of 1300 s. 

The coadded Coma images for each chip were domi- 
nated with light from a number of bright elliptical galax- 
ies. To remove this light, we used elliptical isophotcs 
to construct photometric models of the six brightest 
galaxies on the three WF chips. We then subtracted 
these models from the coadded images. We trimmed the 
elliptical-subtracted coadded images to remove the re- 
gions affected by edge effects, and trimmed the ellipse 
model images to match. We then ring median-filtered 
the elliptical-subtracted and trimmed images to create a 
smoothed map of the background light for each WF chip. 
We subtracted these background light maps to create fi- 
nal background-subtracted images. We also added the 
background light map to the trimmed ellipse model im- 
age for each chip to create final models of the background 
for each chip. 

2.2. Control Field Data 

For control field data, we obtained F606W and 
F814W images of the Hubble Deep Field North (HDF) 
(Williams et al.lll99l from the STScI website 3 . These 
images had been registered, drizzled onto an image with 
a sampling of 0.04" pix -1 , background subtracted and 
normalized to an exposure time of 1 s. 

In the method of differential counts, it is imperative 
to recreate the detection characteristics of the data field 
as closely as possible in th e control field ( for a good 
discussion of this topic , see iBernstein et all l)1995f) and 
lAndreon fc Cuillandrd l)2002j) '). We imposed the Coma 
pixel scale on the HDF frames by rebinning the images 
by a factor of 2.5. To match the Coma images' effective 
exposure time of 1300 s, we multiplied each HDF image 
by 1300. We also trimmed the HDF images to match the 
area of the trimmed Coma images. 

To match the noise of the Coma images, we generated 
a spatially- dependent noise frame for each chip. At each 

3 http: //www. stsci . edu/f tp/science/hdf /hdf .html 



pixel of the noise frame, the expected noise level of the 
corresponding Coma background model pixel was calcu- 
lated based on the background level at that pixel, the 
read noise, and the number of frames that were coadded 
to create the image. A value for the noise pixel was then 
drawn at random from a Gaussian distribution centered 
on zero with standard deviation equal to that expected 
noise value. The generated value was multiplied by 0.85 
a priori to better match the measured noise levels in 
the background-subtracted Coma images. We added the 
resulting noise frames directly to the corresponding re- 
binned, scaled and trimmed HDF chips, as the existing 
Poisson noise in the HDF images is negligible. 

Finally, we subtracted the background of the HDF im- 
ages in the same manner as the Coma images to keep all 
processing steps as similar as possible. 

3. CREATING THE OBJECT CATALOGS 

3.1. Globular Cluster Contamination 

One of the advantages of using HST data is its su- 
perior resolution over ground-based data. This is impor- 
tant in this work because of the danger of globular cluster 
blends contaminating the LF. Globular clusters are point 
sources at the distance of Coma, but two or more globu- 
lar clusters (or other stellar objects) can appear to be a 
single extended object when close enough together, and 
thus be miscl assified as a cluster gal a xy. This effect was 
discussed bv lAndreon fc Cuillandrd (|2002T> . who noted 
that a large number of faint, extended objects in their 
CFHT images of the Coma cluster were resolved into 
separate point sources in corresponding HST images. 

Examination of our preliminary catalogs revealed that 
while globular clusters were certainly visible as separate 
point sources, these separate sources would still some- 
times be blended into a single object by our detection 
algorithm. We found this to be a function of the choice 
of convolution kernel; when the image was convolved 
with a moderate-sized Gaussian kernel (FWHM ~ 3 pix) 
to enhance detection of extended sources, nearby point 
sources could be blended into a single extended object 
(Figure Convolving with a smaller kernel (FWHM 
~ 1.5 pix) prevented this effect, but using such a small 
kernel caused some of the more extended faint galaxies 
to be missed. 

To prevent globular cluster blends from contaminating 
our catalog while still including as many cluster galaxies 
as possible, we chose to employ a two-pass detection and 
photometry procedure. We convolved the images once 
with a small kernel to detect every point source as a 
separate object. We separated the resulting catalog into 
"stars" and galaxies on the basis of central concentration, 
and created a stellar mask for each chip. We applied 
this mask in a second pass with a larger kernel, such 
that pixels associated with stellar objects in the first pass 
could not be detected on the second pass. The results of 
this second pass then formed our catalog. 

The HDF has many fewer point sources than the Coma 
images, and therefore is much less affected by the dan- 
ger of point source blends. However, in the interest of 
processing the data and control images as identically as 
possible, we applied the same two-pass procedure to the 
HDF images as well. 

3.2. Detection and Photometry 
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TABLE 1 

Previous studies of Coma's luminosity function 



Study 


Passband 


Field size 
(arcmin 2 ) 


Pixel sc'tilt' 
f fircsGC pix J 


Lim. rcicLg. 

l V ^-.licL lllcl^io j 




Thompson & Gregory 1993 


b 


14 292 


18.56 b 


20 


-1.43 


Biviano et al. 1995 


b 


2 496 


67.2 b 


20 


-1.2 ± 0.2 


Bernstein et al. 1995 


n 


OZ.A 




ZO.O 




Seeker et al. 1997 


R 


~ 700 


0.53 


22.5 


-1.41 ± 0.05 


Lobo et al. 1997 


V 


1 500 


0.3145 


21 


-1.8 ± 0.05 


Trentham 1998b 


R 


674 


0.22 


23.83 


-1.7 


Adami et al. 2000 


R 


52.2 


0.473 


22.5 


-1 


Andreon & Cuillandre 2002 


B 


720 


0.206 


22.5 


-1.25 




V 


1 044 


0.206 


23.75 


-1.4 




R 


1 044 


0.206 


23.25 


-1.4 


Bcijcrsbergen et al. 2002 


U 


4 680 


0.333 


21.73 


i o 9 +0.018 
± - oz -0.028 




B 


18 720 


0.333 


21.73 


-, o 7 +0.024 

^'-o.oie 




r 


18 720 


0.333 


21.73 


-1 16+ ' 012 


Mobasher et al. 2003 


R 


3 600 


0.21 


19.5 


i io+0.04 
-J--1S-0.02 


Iglesias-Paramo et al. 2003 


r' 


3 600 


0.333 


20.5 


_1 47+0.08 



a When no error is given for a, no value was supplied in the text* 5 These studies were done using photographic 
plates. The value listed for "Pixel scale" is actually that for plate scale, measured in arcsec mm - lc Bernstein 
et al. reported two limiting magnitudes for their work. The first, shown here, is the faintest magnitude bin 
included in the fit to their luminosity function. Beyond this, they believed their counts were accurate to R 
~ 25.5, but decided not to use the data for fear of globular cluster contamination. 




functi 



n 



F(< r) 



(1) 



Fig. 1. — Blended objects in the Coma catalog. The pairs of 
point sources circled above are easily resolvable by eye, but were 
blended into a single extended object by our detection algorithm 



We used the SExtractor package iJBertin fc Arnoutsl 
1996) for detection and photometry. We performed the 
first detection/photometry pass on the F606W images, 
convolving with a FWHM = 1.5 pix Gaussian kernel. We 
set the detection threshold to 1.48 in units of background 
standard deviation in the unconvolved image; this corre- 
sponds to a threshold of ~ 3.5 a on the convolved image. 
We gave the final Coma background model to SExtrac- 
tor as a variance-type weight map for both the Coma 
image and the corresponding HDF image. SExtractor 
uses these weight maps, which describe the noise inten- 
sity at each pixel in units of relative weight or relative 
variance, to slightly adjust the threshold at each pixel to 
compensate for noise varying across the frame. 

To identify the stellar objects in the detection lists, 
we calculated a "simplified" Petrosian radius r petros for 
each object. This is defined as the radius at which the 



reaches its maximum, where F(< r) is the flux within 
an aperture with radius r. Note that rf is proportional 
to signal-to-noise when an image is sky -noise dom i nated . 
This measure is loosely based on the iPetrosianl l)1976fl 
radius, which has been used to determine central con- 
centration for various purposes in numerous studies (see, 
for example, Strauss et al. 2002). 

We plotted the r pe t ros for each object versus its 

MAG BEST magnitude (Figure EJ). We chose a cut of 

fpetros = 1-35 to separate stellar objects from galaxies. 
Some stars do fall above this line; however, as the ob- 
ject of this step is to block as many stars as possible 
while not losing any galaxies, a lower value is preferable. 
We then used the segmentation image produced by SEx- 
tractor to identify the pixels associated with each stellar 
object. We created a weight map where these pixels were 
assigned a weight of zero, while the remaining pixels re- 
ceived a weight equal to the inverse of the value of the 
corresponding Coma background model pixel. 

We performed the second detection/photometry pass 
using this weight map. The image was convolved with a 
FWHM = 3 pix Gaussian kernel, and we set the detection 
threshold to 0.86 (again equivalent to ~ 3.5 sigma on 
the convolved image). We set the weight threshold such 
that no pixels with weights of zero would be included 
in the detected objects. Photometry was performed on 
both the F606W and F814W images; object positions 
and image moments were determined on the (deeper) 
F606W frames, and then applied to both the F606W and 
F814W frames to calculate magnitudes. Fixed-aperture 
magnitudes were also determined in both filters for later 
use in calculating colors. 

3.3. Galaxy Selection 

We trimmed the preliminary Coma catalogs to remove 
objects which fell in the cores of the subtracted Coma el- 
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Instrumental F606W magnitude Instrumental F606W magnitude 



Fig. 2. — Simplified Petrosian radius versus instrumental F606W 
magnitude for objects detected by convolving with a Gaussian ker- 
nel of FWHM = 1.5 pix. For the purposes of constructing the 
globular cluster mask, objects with r pe tros < 1-35 are classified as 
stellar. Some stars do fall above this line; however, a lower value 
is preferable at this step to block as many stars as possible while 
not losing any galaxies. 

liptical galaxies. The HDF had no large elliptical galax- 
ies to subtract; however, its noise was generated from 
the Coma background models containing these galaxies, 
and so we removed HDF objects in the trim regions from 
the catalogs as well. Finally, to remove any remaining 
stellar objects from the catalogs, we cut all objects with 

Rpetros — E5d. 

3.4. Limiting Magnitude 

We performed add-galaxy experiments to determine 
the limiting magnitude of our catalogs. For each of 12 
0.5 magnitude bins, we generated 100 fake galaxies and 
added them to one chip of the Coma and HDF images. 
After running SExtractor on the added-object images, 
we determined the number of added objects that were 
recovered. We repeated this process 20 times for both 
Coma and the HDF and combined the results. 

Figure[3]shows the results for both Coma and the HDF. 
As can be seen, the Coma catalog is over 80% complete 
to an instrumental F606W magnitude of -4.25. This can 
then be considered the limiting instrumental magnitude 
of the catalog. 

3.5. Magnitude Transformation 

Finally, we transformed our magnitudes to the Vega 
system Johnson- Cousins B, R and I bands. To trans- 
form from instrumental magnitudes to the Vega system, 
we used the HDF photometric zeropoints on the STScI 
webpage 4 ; to transform from F606W to B and R, and 
from F814W to I, we used the data of iFukueita et all 
l)1995|) . We determined the F606W - F814W color nec- 
essary for this transformation from the fixed-aperture 
magnitudes calculated by SExtractor. To those objects 
for which SExtractor could not determine fixed-aperture 
magnitudes, we assigned a color of F606W - F814W = 
0.8, chosen as the typical color from a histogram of ob- 
ject colors in the final Coma catalog. Under this trans- 
formation, our limiting magnitude (80% completeness) 



Fig. 3. — Average percentage of added objects recovered in each 
magnitude bin. Open circles represent the results of objects added 
to the Coma field; filled squares are the results of objects added to 
the HDF field. Each point represents the average of 20 runs. Error 
bars are the standard deviation in the mean. 

becomes R = 25.75 in the Vega system. 

3.6. Color and Morphology 

The method of statistical background subtraction can 
be improved when combined with other forms of back- 
ground subtraction. A common technique is to pre-select 
cluster members based on morphology or color before 
performing the statistical background subtraction. This 
reduces the number of background galaxies that must be 
removed statistically, therefore reducing the errors inher- 
ent in statistical subtraction. 

Unfortunately, the objects found in our Coma fields are 
too small to derive meaningful morphology information. 
The objects in the final Coma catalog have a typical ef- 
fective radius of ~ 2 pix. This is simply too few pixels 
to determine morphological type. 

Pre-selecting based on color also proved impossible 
with these data. Figure 0] shows color-magnitude dia- 
grams of objects in the final Coma and HDF catalogs. 
As can be seen, the distributions are very similar. There 
is no obvious way to select Coma cluster members based 
on their color. 

3.7. Cosmic Variance 

Cosmic variance - the fact that background counts vary 
from pointing to pointing due to clustering and large 
scale structure - must be considered when using statis- 
tical background subtraction. To understand the effects 
of cosmic variance on the luminosity function, we evalu- 
ated the expression of lPeeblesI 1)1980(1 for the variance of 
the count N of objects in a randomly placed cell. This 
expression is given by 

((N - 7]fl) 2 ) = rjfl + r] 2 J dOid0 2 w(^i 2 ) (2) 

where 77 is the mean density of objects on the sky, f2 is 
the size of the cell, and w(0 12) is the two-point angular 
correlation function. 

4 http: //www . stsci . edu/f tp/science/hdf /logs/zeropoints . txt 
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Fig. 4. — Color-magnitude diagrams of the final Coma and HDF 
catalogs. Filled circles represent objects from the Coma catalog; 
crosses represent objects from the HDF catalog. 

The two-point angular correlation function of galax- 
ies is usually parameterized by w{9\2) = A w 9~ s , where 
A w is often a function of magnitude. A search of the 
literature yielded no applicable studies of the two-point 
galaxy angular correlation function in th e Vega R band . 
We chose therefore to use the results of I Wilson! (2003), 
who did her work in the Vega V band. Wilson's data 
can be applied to this work if a constant Vega V-R color 
is assumed for the objects in the Coma and HDF cat- 
alogs. As discussed earlier, the typical Vega F606W- 
F814W color of objects in th e final Coma catalog i s 0.8. 
Using the transformations of lFukugfta et all l)1995[) . this 
corresponds to a typical Vega V-R color of 0.5. There- 
fore, Wilson's coefficients for the bin mi < V < m 2 were 
applied to the bin mi — 0.5 < R < mi — 0.5 in this work. 
The values used can be found in Tabled 

For r), we determined the mean density of background 
galaxies from the compilation of published field number 
counts available from the Durham Cosmology Group 5 . 
Because counts in the Vega V band were not available, 
the counts in the Vega R band were used. Following the 
reasoning used above, the counts in the bin mi — 0.5 < 
R < m 2 — 0.5 were used as the counts in the bin mj < 
V < m 2 ■ From the given data, the average count in each 
0.5 mag bin was determined. This value was divided by 
3600 to get the average number of galaxies per arcmin 2 
in each bin. As the cosmic variance was to be calculated 
in 1.0 mag bins, the results from pairs of bins were added 
together. The resulting values can be found in Table |3 

Finally, for fi, we simply calculated the total field size 
of this work by adding together the trimmed areas of the 
three WF chips. This value can also be found in Tabled 

Using these values, we then evaluated Peebles' integral 
numerically for each of the four 1.0 magnitude bins. The 
distance between points was calculated as the linear dis- 
tance; this approximation is valid for fields as small as 
the one used here. The resulting values of cosmic vari- 
ance and standard deviation for each of the bins can be 
found in Table 

5 http://star-www.dur.ac.uk/ nm/pubhtml/counts/counts . html 



To convert the cosmic variance to an error in the lumi- 
nosity function, two steps were necessary. First, as the 
luminosity function was to be expressed in terms of 0.5 
mag bins, the error in each 1.0 mag bin due to cosmic 
variance had to be divided between two bins. We chose 
to do this based on the relative population of the two 
bins: If <tt is the standard deviation due to cosmic vari- 
ance for a 1.0 mag bin, and ni and n 2 are the background 
counts in two 0.5 mag bins, then the standard deviation 
due to cosmic variance in the two 0.5 mag bins is given 

by 

Ol=C7 T ; (3) 

ni + n 2 

a 2 =a T ; 4) 

ni + n 2 

Second, the error due to cosmic variance only applies 
to background counts. For the HDF, all the counts in 
each bin are due to background galaxies. For Coma, 
however, only a fraction of the counts in each bin are 
from background galaxies. The exact number for each 
bin is unknown; statistically, it is assumed to be the same 
as the HDF count for the corresponding bin. Therefore, 
we set the error due to cosmic variance for each Coma 
bin equal to the error due to cosmic variance for the 
corresponding HDF bin. 

We are not able to give an estimate of the cosmic vari- 
ance all the way down to our limiting magnitude, because 
Wilson's study does not go this deep. Therefore, we are 
forced to assume that the error due to cosmic variance 
remains small. Further, by using Wilson's work, we are 
assuming that the conditions of her study - instrument 
sensitivity, seeing, detection and photometry algorithm, 
and so on - are sufficiently similar to the conditions of 
this study to warrant the comparison. 

Obviously, these are very strong assumptions. Ideally, 
we would evaluate cosmic variance by obtaining addi- 
tional data taken and reduced under the same conditions 
of the Coma data. However, this is beyond the scope of 
this paper, and so we continue under the assumptions 
stated above. 

4. RESULTS 
4.1. Comparison to Other Dwarfs 

To compare the objects in our Coma catalog to known 
dwarf galaxies, we collected six published catalogs of 
dwarf galaxies (Table[3J| and obtained the magnitude and 
half-light radius for each of the galaxies in physical units. 
When a distance to the group or galaxy was given with 
the catalog, we adopted that value. For the Virgo cluster, 
we followed iKavelaars et al.l 1(2000]) and used a distance 
modulus of (m-M) = 30.99. 

We converted the magnitudes and half-light radii 
of objects in our final Coma catalog into physical 
units, adopting a dist ance modulus of (m-M) = 35.05 
l|Kavelaars et al.ll2000(l . For comparison, we also deter- 
mined magnitudes and half-light radii for Coma's glob- 
ular clusters (i.e., the objects found to be stellar while 
constructing the globular cluster mask). 

Figure [S] shows half-light radius as a function of mag- 
nitude for all the objects. Although the Coma galaxies 
are much fainter than most known dwarf galaxies, they 
fall in the same general locus. The distribution is under- 
standably more scattered: background galaxies are also 
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TABLE 2 
Cosmic Variance 



Vega 


R bin 


Vega 


V bin 


5 


logio^MT) 


V 


n 


Cosmic 


Standard 














(arcmin~ 2 mag -1 ) 


(arcmin 2 ) 


Variance 


Deviation 


20.5 


- 21.5 


21.0 


- 22.0 


0.8 


-1.13 ± 0.06 


1.15 


4.62 


7.72 


2.78 


21.5 


- 22.5 


22.0 


- 23.0 


0.8 


-1.49 ± 0.05 


2.69 


4.62 


21.3 


4.62 


22.5 


- 23.5 


23.0 


- 24.0 


0.8 


-1.72 ± 0.05 


6.37 


4.62 


48.4 


6.96 


23.5 


- 24.5 


24.0 


- 25.0 


0.8 


-2.26 ± 0.09 


14.6 


4.62 


96.2 


9.81 



TABLE 3 

Published catalogs of dwarf galaxies 



TABLE 4 
The luminosity function 



Catalog 



Number of galaxies 



Reference 



m R 



N 



AN P ANcv log a Alog cr 



Local Group 
M101 
M81 
Northern Field 
Southern Field 
Virgo Cluster 



21 
18 
19 
17 
21 
365 



Bender et al. 
Bremnes et al. 
Bremnes et al. 
Barazza et al. 

Parodi et al. 



1992 
1999 
1998 
2001 

2002 



Binggeli & Cameron 1993 



present in the Coma catalog, and these show a range of 
half-light radii at each magnitude. The Coma globular 
clusters, however, clearly occupy a very different region 
in the plot. From this plot, we conclude that the ob- 
jects in the Coma catalog are not globular clusters, and 
do have the same general morphology as known dwarf 
galaxies. 

4.2. The Luminosity Function 

Table0]shows the data that were used to construct the 
LF. The number of galaxies in the HDF image is scaled to 
the unblocked area of the Coma image; the Coma image 
has many more point sources, and therefore more pixels 
are blocked by the stellar mask. 

Poisson counting errors and errors due to cosmic vari- 
ance are given separately. The error due to cosmic vari- 
ance in each 0.5 mag bin could be added in quadrature to 
the Poisson error. However, the calculation here was per- 
formed only to serve as an estimate of the effects of cos- 
mic variance on the slope of the luminosity function. A 
full treatment would have to account for the correlations 
in error between magnitude bins introduced by cosmic 
variance. Due to the simplistic approach and many as- 
sumptions used here, we choose to present the two errors 
separately. 

Figure [fj] shows the number counts of the Coma field 
and the scaled HDF field. As can be seen, cluster counts 
do not start to dominate over the background until the 
last two bins. This is due to the small field size surveyed 
here: cluster galaxies with 20 < R < 24 are so rare that 
< 1 is expected in a field this size. The galaxies in this 
magnitude range that are detected in the Coma field are 
due to background contamination, and so their numbers 
match those found in the background-only field. 

Figure [7| shows the result of subtracting the HDF 
counts from the Coma counts - the luminosity function. 
The shape of the luminosity function is generally param- 
eterized by the Schechter function, which can be written 
in terms of magnitude as 

4>(M)dM ~ i -0-4(«+i)A/ e -io°-<«*--» (5) 

At faint magnitudes, the first term dominates and so the 
slope of the logarithmic luminosity function is directly 
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Note. — The columns arc: rriR — midpoint of bin 
in apparent Vega R magnitudes; N — number of galax- 
ies; AN P - Poisson error in the count; ANcv ~ error 
due to cosmic variance in the count; log a — logarith- 
mic surface density of objects (number of galaxies per 
deg 2 ); Alog a — Poisson error in log a. 



related to the parameter a. Fitting a straight line to 
our logarithmic LF using a weighted least-squares fit, we 
obtain a = -2.29 ±0.33. 

5. DISCUSSION 
5.1. Comparison to Other Work 
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Fig. 5. — Logarithm of effective radius in kpc versus absolute Vega B magnitude for various dwarfs. Data for dwarfs in Coma are from 
this work; references for the other catalogs used can be found in Table [3] Also shown for comparison are the Coma globular clusters 
identified in this work. 



We have found an extremely steep slope for the faint 
end of the luminosity function in the Coma cluster. This 
is especially striking when compared to studies of the 
faint-end LF in clusters such as Virgo and Fornax. These 
clusters are close enough that measurements can easily 
be made at the same absolute ma gnitudes reached here. 
For ex ample, in the Virgo cluster. TTrentham fe Ho dgkin 
(2002) found magnitude-dependent slopes ranging from 
-1.49 ± 0.11 < a < -1.02 ± 0.09 for -15 < M B < 

— 12 (roughly equivalent t o our absolute m agnitude range 
of -13 < M R < -9). iSabatini et all 1)200311 found a 
slightly steeper slope of a = —1.6 for Virgo in the range 

— 14 < Mb < —10, though thi s is still significantly flatter 
than our slope. For Fornax, Milker et alJ l)2003|) found 
an extremely fiat slope of a = 1.1 ± 0.1 in the range 
-16 < M v < -9. 

Another interesting comparison is to look at other 
studies of the LF in the Coma cluster. As noted in the 
Introduction, many other such studies have appeared in 
the literature. Tabled indicates that little consensus has 



been reached on the "true" slope of the LF. To better 
compare these studies, we constructed a composite LF. 

We obtained background-subtracted number counts for 
the Coma cluster from the studies listed in Table QJtliat 
were conducted in the R band. We converted the magni- 
tudes to absolute ma gnitudes using a distance modulus 
of (m - M) = 34.83 i|Trenthamlll998b() . We then scaled 
the counts and errors to an area of 1 deg 2 , using the field 
size reported by the authors. 

No other normalization was performed. This intro- 
duces some scatter to the composite LF; to be perfectly 
correct, larger surveys - covering more of the low density 
cluster outskirts and hence having a lower average sur- 
face density of galaxies - should be scaled differently than 
surveys covering only the high density cluster core. How- 
ever, taking the scale radius (i.e., the radius at which the 
slope of the profile is the average of the inner and outer 
slope) of the Co ma cluster to be r* = 32 kpc for H G = 70 
km s _1 Mpc -1 llLokas fc MamorJl2003D. the largest sur- 
vey included here l)Moblish^r^tniLll200 3i) covers an area 
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Fig. 6. — Number counts from the Coma cluster field and the 
HDF field (after scaling). Error bars show the Poisson error in the 
counts. 
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Fig. 7. — Vega R band luminosity function of the Coma cluster, 
expressed as logarithmic counts per 0.5 magnitude bin. Error bars 
show the Poisson error in the counts. The line is a weighted least 
squares fit to the points. 

roughly 9 times that of the Coma cluster core. Assuming 
a standard surface density profile (cr(r) = (l+(r/r s )™) _1 , 
for n = 1 or 2), moving to a radius ~ 3r s would roughly 
halve the average surface density determined from ob- 
serving only the cluster core. This corresponds to a scat- 
ter of ±0.3 in our logarithmic composite LF, an accept- 
able level of error for our purposes. 

Figure [S] shows the composite R band luminosity func- 
tion for the Coma cluster. The bright end shows good 
agreement among the various surveys. At the faint end, 
the counts of this work are a clean continuation of the 
brighter counts from the literature and clearly show a 
steep faint end slope for the Coma clust er LF. 

Esp eciall y interesting a re the counts o flMobasher et all 
( 2003) and lAdami eTaf] l|2000]) . Of the R band studies 
discussed here, they were the only groups to use spec- 
troscopic redshifts to determine cluster membership. As 



can be seen in Tabled their values for a are among the 
lowest published. This trend of lower a values for spec- 
troscopic studies has sometimes been interpreted as an 
indication that statistical background subtraction causes 
artificially steep slopes. 

However, looking at the composite LF, the spectro- 
scopic counts of Mobasher et al. agree very well with 
those obtained by using photometric methods to deter- 
mine cluster membership. Their counts stop before mov- 
ing into the truly faint regimes of the LF; spectroscopic 
studies, by their nature, are limited to bright galaxies 
and so rarely probe the faint end of the LF. The flatter 
slopes reported by spectroscopic studies may simply be 
a result of probing only the brighter, less steeply rising 
part of the LF. 

The counts of Adami et al. do extend into fainter re- 
gions of the luminosity function. Their study worked 
with a small sample of only 88 redshifts. They did not 
calculate error bars for their luminosity function, but as- 
sumed them to b e the same as those show n for the (much 
higher) counts of iBernstein et all l)1995j) . With these er- 
ror bars, their counts are consistent at a ler level with 
those obtained using photometric methods. The general 
trend of the points may hint at a possible systematic er- 
ror arising from using photometric methods to construct 
a luminosity function, but it is dangerous to draw con- 
clusions based on such small numbers and rough error 
estimates. Unless a larger, better characterized sample 
of redshifts are obtained which show the same results, we 
see no compelling evidence that statistical background 
subtraction leads to an artificially steep LF slope. 

Perhaps the most important thing that Figure [5] shows 
is that, although the various studies did not agree on a 
value for a, the counts themselves agree very well. Cau- 
tion must always be used when comparing results via 
parameterized values rather than the underlying data. 

5.2. Implications of a Steep Faint-End Slope 

We have found a steep faint end slope to the Coma 
LF. Such steep slopes are not seen in all environments. 
What could cause this difference between faint end slopes 
in different environments? 

Qualitatively, o ur results agree with the idea of 
iTullv et al] (j2fl02): Coma, being a rich cluster, formed 
many dwarf galaxies before reionization "squelched" 
dwarf galaxy formation. A more quantitat i ve an alysis of 
this idea was performed bv iBenson et al.l l)2003j) . They 
used detailed semi-analytic models of galaxy formation, 
including the effects of supernova feedback, photoion- 
ization suppression, and dynamical friction and mass 
loss due to tidal forces, to determine the properties of 
galaxies in a range of environments. They found that 
when photoionization suppression was switched on, faint 
end slopes in poor environments did become flatter than 
those in rich environments. Howev er, comparing their re- 
sults t o the observational data of iTrentham fc Ho dgkin 
(2002), they found that photoionization alone cannot 
flatten the slope enough to match the luminosity func- 
tions observed in the Local Group and Ursa Major. 

Thi s conclusion is supported by iGrebel k, Gallagher! 
( 2004) , who examined the ages of stellar populations in 
Local Group dwarfs. They found no evidence for a halt 
or reduction in star formation at 20 < z < 6.4, the typi- 
cal range of epochs for reionization. This indicates that 
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Fig. 8. — The composite Vega R band luminosity function of the Coma cluster, expressed as logarithmic counts per 0.5 magnitude bin 
and per deg . Bernstein et al.'s "additional counts " are those not included in their final luminosity function for fear of globular cluster 
contamination. All counts arc scaled to an area of 1 deg 2 ; no other normalization was performed. 



rcionization cannot be the dominant influence in the evo- 
lution of Local Group dwarfs. Other feedback mecha- 
nisms must also be at work. 

Obviously, further theoretical efforts are needed to re- 
fine the exact processes that shape the luminosity func- 
tion. Based on the results of this and similar work, these 
studies must now also explain how these processes can 
lead to a luminosity function that is different in environ- 
ments of different densities. 

6. CONCLUSION 

We have found a steep faint end slope for the galaxy 
luminosity function in the core of the Coma cluster. Us- 
ing the method of statistical background subtraction and 
archival HST images has allowed us to achieve a limit- 
ing magnitude of R = 25.75, making this the faintest 
and highest resolution determination of the Coma lu- 
minosity function to date. The fine pixel scale of the 
WFPC2 camera and our globular cluster mask enabled 
us to effectively eliminate globular cluster contamination 



from the final catalog. We find the slope of the luminos- 
ity function to be fit by a = -2.29 ± 0.33. While this 
value is not in agreement with other published values of 
a for the Coma cluster, a composite luminosity function 
shows our counts to be in good agreement with previous 
counts. This result could be affected if the errors due to 
cosmic variance are larger than our estimate. Assuming 
they are not, the steep slope found in this work is consis- 
tent with theories that predict photoionization and other 
feedback effects will affect environments of low density 
more severely than environments of high density. 
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